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HIGH-LATITUDE Hi IN THE LOW SURFACE BRIGHTNESS GALAXY UGC 7321 



m 
O 
O 

(N 



(N 



> 
(N 
O 

o 
m 
o 

Or 

6 ■ 

in 
c3 



L. D. Matthews 1 
Kenneth Wood 2 
Draft version February 2, 2008 

ABSTRACT 

From the analysis of sensitive Hi 21-cm line observations, we find evidence for vertically extended Hi 
emission (\z\ <2.4 kpc) in the edge-on, low surface brightness spiral galaxy UGC 7321. Three-dimensional 
modelling suggests that the H I disk of UGC 7321 is both warped and flared, but that neither effect can 
fully reproduce the spatial distribution and kinematics of the highest z-height gas. We are able to model 
the high- latitude emission as an additional H I component in the form of a "thick disk" or "halo" with 
a FWHM~3.3 kpc. We find tentative evidence that the vertically extended gas declines in rotational 
velocity as a function of z, although we are unable to completely rule out models with constant V(z). In 
spite of the low star formation rate of UGC 7321, energy from supernovae may be sufficient to sustain this 
high-latitude gas. However, alternative origins for this material, such as slow, sustained infall, cannot 
yet be excluded. 

Subject headings: galaxies: ISM — galaxies: spiral — ISM: structure — galaxies: halos — galaxies: 
individual (UGC 7321) 



1. BACKGROUND 

UGC 7321 is a nearby, isolated Sd spiral galaxy seen 
nearly edge-on to our line of sight (i w 88°). The dis- 
tance to UGC 7321 is a matter of some uncertainty, but 
based on the brightest resolved stars it is estimated to 
be «10 Mpc (see Uson & Matthews 2003), and we adopt 
this value in the present work. UGC 7321 exhibits a 
highly flattened, diffuse stellar disk with no bulge com- 
ponent (Figure 1), and a very small intrinsic stellar scale 
height (h z ~150 pc; Matthews 2000). Its photometric 
properties have been studied by Matthews, Gallagher, 
& van Driel (1999), who found that it has an intrinsi- 
cally low optical surface brightness disk [/ig,j(0) «23.5 
mag arcsec -2 ], indicative of a low current star formation 
rate. This is corroborated by its low far-infrared lumi- 
nosity (Lpm ~ 7 x 10 7 L Q ), the nondetection of 21-cm 
radio continuum emission by Uson & Matthews (2003; 
Fcont < 0.41 ±0.25 mJy), and the rather modest Ha emis- 
sion in this galaxy (Lu a ~ 1 x 10 40 erg s -1 ; see Sec- 
tion 4.2). The high neutral gas content of UGC 7321 rel- 
ative to its blue luminosity (AIhi / A4q/£q), to- 
gether with its low dust content (Matthews & Wood 2001), 
moreover suggest that this galaxy has been an inefficient 
star-former over most of its lifetime. 

Uson & Matthews (2003) recently obtained new aper- 
ture synthesis imaging observations of UGC 7321 in the 
Hi 21-cm line using the Very Large Array (VLA) 3 . De- 
tails regarding these observations, as well as the global and 
kinematic properties of the neutral hydrogen in UGC 7321 
are described by those authors. With a total of 12.3 hours 
of on-source integration time, the VLA data for UGC 7321 
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are among the most sensitive H I observations of edge-on 
spiral galaxies obtained to date. Since the observations 
also have reasonably good spatial and spectral resolution 
(~ 16" and ~5.2 km s _1 , respectively), they are well- 
suited for structural and kinematic studies of the H I gas 
in this nearby system. Here we use these data to explore 
the vertical distribution of the neutral hydrogen gas in 
UGC 7321. With the aid of three-dimensional models, 
we report evidence that the structure and kinematics of 
the Hi component in UGC 7321 are more complex than 
a simple, thin-disk distribution. UGC 7321 has a cool Hi 
disk that appears to be both warped and flared, but our 
analysis also suggests the presence of additional vertically- 
extended emission (Section 3.3). In Section ?? we briefly 
discuss possible origins for this high-latitude material. 

2. HI HALOS OF GALAXIES 

H I surveys of the Milky Way have revealed that the dis- 
tribution of its neutral hydrogen component is far more 
complex than simply a cool, thin disk of material ly- 
ing along the Galactic midplane. The Galactic neutral 
medium instead contains both cool and warm layers of dif- 
fering scale heights (Lockman 1984), as well as a halo com- 
ponent extending to substantial z-heights (h z ^4.4 kpc; 
Kalberla et al. 1998). 

From sensitive aperture synthesis observations, high- 
latitude H I gas has now also been detected in other nearby 
galaxies. For example, using deep Hi observations of 
the nearby, edge-on Sbc spiral NGC 891, Swaters, San- 
cisi, & van der Hulst (1997) reported evidence for Hi gas 
extending to heights of at least 5 kpc above the plane. 
Based on three-dimensional modelling, these authors ar- 
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gued that this emission is part of an H I halo that lags in 
rotation speed by 25-100 km s _1 relative to material in 
the disk. 4 Similar analyses for the moderate inclination 
spiral NGC 2403 led to the suggestion that this galaxy 
also has a rotationally lagging Hi halo and anomalous 
H I material that may be related to Galactic high- velocity 
clouds (Schaap, Sancisi, & Swatcrs 2000; Fraternali et al. 
2001,2002). In addition, high-latitude Hi features now 
have been seen in NGC 5775 and NGC 2613 (Irwin 1994; 
Irwin & Chaves 2003). Interestingly, NGC 5775 shows evi- 
dence for an Hi halo, while NGC 2613 seems to show only 
multiple discrete high-latitude features (Irwin & Chaves 
2003). 

For NGC 891 and NGC 2403, Swaters et al. (1997) and 
Schaap ct al. (2000) , respectively, postulated that the H I 
halo gas in these galaxies is likely to have originated from 
a galactic fountain-type mechanism of the type described 
by Shapiro & Field (1976) and Bregman (1980). This ex- 
planation seems plausible for these systems, since both are 
bright galaxies with moderately high star formation rates 
and significant Ha and radio continuum emission. Still, 
the sample of nearby galaxies for which gaseous halos, 
particularly Hi halos, have been detected remains quite 
small, and a number of puzzles regarding their structure 
and origin remain. 

One difficulty with the galactic fountain scenario for the 
origin of high-latitude gas is that the energetics required 
to raise material from the disk to z-heights of several kpc 
often appears to be much larger than can seemingly be 
supplied by ordinary supernova events alone (e.g., Rand 
& van der Hulst 1993; Lee & Irwin 1997). Moreover, ev- 
idence is accumulating that the observed vertical velocity 
gradients in halos are not entirely consistent with fountain 
predictions (Collins, Benjamin, & Rand 2002). 

In some cases, it has been suggested that extraplanar 
gas, including Hi, may have been acquired during interac- 
tions (e.g., Sancisi 1988) or buoyed up by tidally-triggered 
instabilities. Nonetheless, interactions seem an unlikely 
source for the halo gas in cither NGC 891 or NGC 2403, 
since both of these galaxies appear to be isolated. One 
possible alternative for the origin of at least some fraction 
of the high-z H I material in galaxies like these is that this 
gas is a remnant from an earlier evolutionary phase or is 
accreted via a slow "drizzle" (e.g., Oort 1970). The in- 
terest in these types of scenarios is further fueled by the 
recent results from QSO absorption line experiments to- 
ward moderate redshift galaxies (redshifts <1) that have 
shown these systems frequently are surrounded by very ex- 
tended gaseous halos with covering factors of near unity 
(e.g., Lanzetta et al. 1995; Chen et al. 2001b). The 
large spatial extents of these halos (~180 h~ x kpc), their 
presence around galaxies spanning all morphological types, 
and the inferred correlations between halo size and galaxy 
mass (Chen et al. 2001b), all hint that their origin may 
be tied to processes intrinsic to the galaxy formation pro- 
cess rather than the details of current star formation (see 
also Section ?? below). It follows that at least some rem- 
nant of such halos still may be present around many disk 
galaxies in the local universe, particularly if the halos are 
continually replenished. 

One important way to test whether at least some frac- 

4 Benjamin 2002 points out that the high end of this range appears 



tion of the high-latitude H I emission observed in spiral 
disks could be linked to the galaxy formation process 
rather than a by-product of vigorous star formation, tidal 
interactions, or mergers is to search for evidence of high 
z-height Hi emission in isolated disk galaxies with very 
low levels of current star formation. H i-rich, low surface 
brightness (LSB) "superthin" spirals are just such objects 
(e.g., Uson & Matthews 2003). Since high z-height gas can 
be most unambiguously detected in highly-inclined galax- 
ies, we have chosen the edge-on LSB spiral UGC 7321 for 
such an investigation. 

3. RESULTS 
3.1. The HI Total Intensity Map 

A total Hi intensity map of UGC 7321 is shown in 
Figure 1, overlayed on an optical i?-band image from 
Matthews et al. (1999). It is apparent from this figure 
that even after accounting for the finite resolution of the 
VLA beam, the Hi emission in UGC 7321 is not exclu- 
sively confined to a very thin layer. If the Hi disk had 
an intrinsic FWHM thickness less than or comparable to 
that of the stars (i.e., <5"3; Matthews 2000), the ap- 
parent thickness as observed with a 16" FWHM Gaussian 
beam should be ~17". Instead we find based on intensity 
slices parallel to the minor axis that the Hi layer has a 
mean apparent thickness of ~ 24" FWHM. In addition, 
a number of filamentary features protrude to even higher 
z. That these features are not simply noise is implied by 
the fact that the disturbances are evident in several suc- 
cessive underlying contours (see also Figure 7 of Uson & 
Matthews 2003). 

One-dimensional H I intensity slices extracted perpen- 
dicular to the disk of UGC 7321 at various projected radii 
consistently show cores that are Gaussian in shape, but 
with additional extended wings (Figure 2). These wings 
stretch to z-distances of up to ~ ±50" (^2.4 kpc). Quali- 
tatively the profiles shown in Figure 2 resemble plots of 
the distribution of H I volume density as a function of 
z in the Milky Way disk (see Lockman 1984, his Fig- 
ure 5). In the case of the Milky Way, Lockman showed 
that the observed density profile could be modelled as the 
sum of two Gaussian components, together with a more 
diffuse and extended exponential component. After cor- 
recting for underestimated low-z emission and a distance 
to the Galactic Center of i?o=8.5 kpc, Dickey & Lock- 
man (1990) give FWHM for the Gaussian components of 
212 pc and 530 pc, and a scale height for the exponential 
component of 403 pc. The nature of the vertical inten- 
sity cuts in UGC 7321 therefore hints at the possibility of 
an analogous multi-component Hi disk structure in this 
"superthin" LSB galaxy. 

While suggestive of an extended gas component, the ver- 
tical intensity cuts and total intensity map that we have 
derived for UGC 7321 give us only a projected density 
distribution and do not permit us to assess at what galac- 
tocentric radius the vertically extended material is located. 
For example, it is feasible that warped or flared gas from 
the outer disk regions seen along the line of sight is af- 
fecting the shape of the wings of the observed intensity 
profiles, and thus could mimic the signatures of a multi- 
component disk. More firmly establishing the presence of 

be inconsistent with Ha measurements of this galaxy. 
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a multiple-component Hi disk in UGC 7321 therefore re- 
quires exploiting the full three-dimensional nature of the 
H I data cube, and utilizing velocity as well as spatial in- 
formation. We undertake this exercise in the following 
subsections. 

3.2. Minor-Axis Position- Velocity Profiles 

One powerful way of probing the structure of verti- 
cally extended gas in disk galaxies is to examine position- 
velocity (P-V) profiles perpendicular to the disk. In Fig- 
ure 3 (right column) we plot such diagrams for UGC 7321 
at several disk locations. In these panels, H I emission can 
be seen extending up to ~ 50" (~2.4 kpc) from the plane. 
At these and other disk locations, most, although not all, 
of this vertically extended material has velocities within 
<50 km s _1 of the systemic velocity. 

In a galaxy disk viewed edge-on, gas that is extended 
in z and that is primarily clustered near the systemic ve- 
locity can be a signature of flaring of the Hi layer (San- 
cisi & Allen 1979). However, certain characteristics of the 
high-latitude material in our present data argue against 
this explanation for UGC 7321. One is the existence of 
asymmetries in these features on the ±z versus — z sides 
of the disk (see Figure 2 & 3), as well on the east and 
west sides (compare the panels at r = ±20" and ±40" in 
Figure 2). For such a highly inclined galaxy, a flare should 
be expected to be largely symmetric about the midplane. 
East-west asymmetries do not strictly exclude a flare, but 
would imply that the midplane density is different on the 
two sides of the disk. A second argument against a pure 
flare model is that the degree of disk flaring required to ex- 
plain the high-z emission would produce unique signatures 
in both the Hi total intensity map and in the individual 
channel maps. As we show below, these features are not 
seen. Thus as we demonstrate in the next sections, while 
UGC 7321 does exhibit manifestations of both warping 
and flaring in its H I disk, neither of these effects alone ap- 
pears to be sufficient to consistently explain the morphol- 
ogy and kinematics of its highest-latitude H I emission. 

3.3. Three- Dimensional Galaxy Modelling 

In order to further probe the nature of the high-z-hcight 
Hi emission in UGC 7321, we have computed full three- 
dimensional (3-D) models of the galaxy and compared 
these with the observed Hi channel maps, vertical P-V 
profiles, and total H I intensity map. For this analysis we 
have modelled the galaxy velocity field using the rotation 
curve for UGC 7321 derived by Uson & Matthews (2003), 
with the approaching and receding sides of the disk aver- 
aged. Assuming a constant Gaussian velocity dispersion 
for the gas of 7 km s _1 , we then computed the Hi inten- 
sity and velocity distribution using a modified version of 
the radiative transfer and rotation curve modelling code 
written by one of us (KW; see Matthews & Wood 2001). 
These model cubes could be subsequently analyzed in a 
manner identical to the real H I data set. 

For our present analysis we have assumed that the Hi 
is optically thin and that its density profile is Gaussian 
in both the r and z directions. The assumption that the 
gas is largely optically thin appears to be justified, since 
Uson & Matthews (2003) found evidence of only very mild 

5 The full set of channel maps for UGC 7321 is presented in Uson & 



self- absorption toward the center of UGC 7321 (where gas 
densities should be highest, and velocity crowding most 
extreme). Moreover, any optically thick clumps are likely 
to be confined primarily to the midplane, whereas in the 
present paper we are mainly concerned with probing the 
spatial and velocity structure of the higher latitude gas. 

The choice of a Gaussian for the vertical intensity distri- 
bution is motivated by its success in describing the bulk of 
the global, vertical Hi distribution of other galaxies [e.g., 
the Milky Way; Lockman (1984)] and its success in fitting 
our present data. However, a sech 2 z distribution (as ex- 
pected for an isothermal disk) has a very similar shape to 
a Gaussian and would also be consistent with the data. 
This choice does not affect our conclusions. 

For UGC 7321 we also find a Gaussian to provide a good 
representation of the radial Hi intensity distribution, as 
is commonly found in late-type spirals (e.g., Krumm & 
Salpeter 1980). Other forms for the radial distribution 
(such as an exponential) provided significantly less satis- 
factory fits. In all of our models, the H I intensity in the r 
direction has a scale length of 5.6 kpc (13.2 kpc FWHM) 
and is truncated at a radius of r mox =11.8 kpc. All of our 
model data cubes also have been convolved with a Gaus- 
sian in the spatial domain in order to match the resolution 
of the VLA data. We adopt an inclination for UGC 7321 
of i = 88° and a position angle of 82°. All of our model 
parameters are quoted as their values prior to convolution 
with the telescope beam. For each of the classes of models 
we describe below, we have explored a wide range of rele- 
vant parameter space, but limit our current discussion to 
only the best-fitting model of each category. 

3.3.1. Single- Component Disk Models 

The first class of model that we considered for UGC 7321 
was that of a smooth, azimuthally symmetric Gaussian in- 
tensity distribution in both the r and z directions (here- 
after, the "smooth model"). We already know this model 
to be somewhat of an oversimplification, since UGC 7321's 
gas disk is clearly is mildly warped (Figure 1; see also Uson 
& Matthews 2003). However, it serves as an important 
baseline for comparison with more complex models. 

For the smooth model, we found the best overall fit to 
the data using a FWHM Hi thickness of 370 pc, and a 
central density of ~0.19 cm~ 3 . In Figure 4 we plot several 
of the channel maps produced from this model alongside 
the real data 5 . In Figure 3 we plot several representative 
minor axis P-V plots for this model. 

We see from Figure 4 that the distribution of the emis- 
sion in most of the smooth model channel maps is not 
a good match to the data; in some of the channels, the 
model emission narrows in z-extent with increasing r (e.g., 
the —72.3 and —98.2 km s _1 panels), while in the actual 
data, the z-distribution of emission broadens with increas- 
ing r. In addition, the smooth model cannot simultane- 
ously match the intensity in the midplane while still re- 
producing the full observed z-extent of the emission (see 
also Figure 3). Finally, the mean FWHM thickness of the 
gas at r = inferred from this model is larger than that 
of the stars in the near-infrared [(FHWM) 4i jv;b «260 pc 
at r — 0; Matthews (2000)], and exceeds typical Hi scale 
heights of normal spiral galaxies. For example, the mean 

Matthews 2003. 
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FWHM of the Milky Way's cool H I disk at intermedi- 
ate galactocentric radii (4-8 kpc) is only ^230 pc (Dickey 
& Lockman 1990). Since in general, all but the youngest 
stars are dynamically hotter than the cool H I disks of nor- 
mal galaxies, this smooth model is unappealing on physical 
grounds, particularly for such a dynamically cold system 
as UGC 7321. Some of these problems can be partially 
alleviated by modelling the emission in terms of two disk 
components of differing scale heights. However, such a 
model still cannot reproduce the broadening of the gas 
in the z direction with increasing projected radius that is 
seen in a number of the individual channel maps. 

Using our smooth model, we have explored the addi- 
tional possibility that this model produces a poor match to 
the data simply because we have overestimated the inclina- 
tion of UGC 7321. However, computing a similar model to 
that described above, except with i — 86°, yields charac- 
teristic "V-shaped" channel maps, unlike those seen in our 
data (see also Figure 4 of Swaters et al. 1997). Moreover, 
this i — 86° model predicts too little density concentra- 
tion along the midplane. We conclude that our adopted 
inclination for UGC 7321 is not significantly in error. 

It is clear from the Hi total intensity map in Figure 1 
that the H I disk of UGC 7321 has a mild integral sign warp 
(see also Uson & Matthews 2003). Therefore a more ac- 
curate model of UGC 7321 naturally should take this into 
account. In addition, inspection of the individual chan- 
nel maps of UGC 7321 suggests that the Hi layer appears 
to be flared (i.e., the z-extent of the gas increases with 
increasing projected radius). In our next set of computa- 
tions, we therefore introduced disk warping and Baring as 
follows: 

p HI oc cxp[-0.5(r/h r ) 2 ] exp[-0.5{z'/h f ) 2 } . (1) 

The flaring scaleheight is hf = h\ + hfor 5 ^ 4 . The warp is 
introduced by setting z' — \z — z warp |, z warp — z w cos<j> 
with (j) the azimuthal angle in the disk midplane, and 
z w = z w0 {r - r w ) 5 / 4 for r > r w . 

This analytic form for the flare was adopted based on 
its ability to provide a good fit to the data. Models where 
the scale height remained nearly constant within the inner 
few kpc (as is inferred, for example, in the Milky Way; 
Burton 1988) produced noticeably poorer fits. 

The radius where the warp commences was taken to 
be 5.8 kpc (the edge of the stellar disk), and its max- 
imum amplitude to be 480 pc (i.e., r w — 5.8 kpc and 
z w o — 0.05 kpc). We assume that the line of nodes of the 
warp is straight throughout the disk. For the flaring, we 
found our best overall fit to the data taking the FWHM 
of the disk to be 260 pc at its center and taking the max- 
imum thickness at r max to be 1.64 kpc FWHM. The cor- 
responding scaleheight parameters are hi = 0.11 kpc and 
h f0 = 0.027 kpc. 

Although it is difficult to accurately constrain gas scale 
heights at small galactocentric radius for edge-on disks 
(e.g., Oiling 1997a,b), we found that models with much 
smaller gas scale heights at r=0 provided less acceptable 
fits to the data. This would imply that the scale height of 
the cool Hi disk in UGC 7321 is not significantly smaller 
than that of the stars. In fact, we note that with our 
adopted parameters, in order to properly fit the outer Hi 
disk scale height, we find that the resulting FWHM of the 



gas layer in our best model actually exceeds by ~ 20% that 
of the stars as measured in the NIR near the outskirts of 
the stellar disk (Matthews 2000). The requirement that 
hgas < therefore provides some estimate of the un- 
certainty in our flare model parameters. Interestingly, a 
near-equality of the gaseous and stellar scale heights in the 
outer disk of UGC 7321 was also suggested by Gallagher 
& Matthews (2002) based on the vertical distribution of 
luminous outer disk stars and dark nebulae seen in Hubble 
Space Telescope images. 

Our resulting warp+flare models are illustrated in Fig- 
ure 3 and Figure 4. Figures 3 & 4 show that our 
warp+flare model results is a substantial improvement in 
fitting the channel maps and the minor axis P-V plots. 
Overall this model reproduces quite well the bulk of the 
brightest midplane emission in UGC 7321. However, we 
find that the warp+flare model still is unable to repro- 
duce the more diffuse, highest z-hcight emission. This 
shortcoming is most evident for channels corresponding to 
velocities within ~60 km s _1 of systemic. Furthermore, 
one-dimensional intensity cuts extracted in a manner anal- 
ogous to those in Figure 2 do not exhibit the same type of 
extended wings seen in the data; rather for the warp+flare 
model, the intensity distribution can be almost completely 
reproduced by a single (albeit broader) Gaussian (Fig- 
ure 5). 

Constraining the peak amplitude of the flare to match 
the maximum observed z-extents of the gas (i.e., FWHM 
of the flare equal to ^3.3 kpc at r max ) significantly lessens 
the agreement with the data, as the highest z-height flared 
material does not follow either the spatial distribution or 
the velocity structure of the observed high latitude gas 
(viz. the "large flare model" in Figures 3 & 4; note partic- 
ularly the —72.3 and —98.2 km s _1 channel maps). More- 
over, models with a higher amplitude flare fail to match 
the observed total H I intensity distribution since the disk 
begins to appear thickest at its edges, contrary to what is 
observed (Figure 6). 

3.3.2. Two-Component Disk Models 

In our own Galaxy, current evidence suggests that 
within the solar radius, Hi material with |z| <1 kpc (or 
possibly as high as |z| <3 kpc) corotates with the gas in 
the plane (Dickey & Lockman 1990; Benjamin 2000; Lock- 
man 2002). However, alternate tracers of the gaseous halo 
at higher z-heights (e.g., Civ) have shown some evidence 
for lagging rotation within the inner Galaxy (Savage & 
Massa 1987). In external galaxies, evidence of lagging ha- 
los of ionized gas has also been reported (Rand 1997,2000; 
Tiillmann et al. 2000). In addition, as noted above, evi- 
dence for rotationally lagging H I halos has now been found 
for a few other carefully studied spirals. We therefore ex- 
plore the possibility that such a model might also describe 
the high-latitude Hi emission structure in UGC 7321. 

For these "lagging halo" models we have included the 
following two components: (1) a thin disk of material that 
rotates according to the disk rotation curve and whose pa- 
rameters (except for the central density) are fixed to those 
derived for our best-fitting warp+flare model above; (2) 
a thicker, Gaussian-shaped layer whose rotational velocity 
is permitted to decrease linearly as a function of z-height. 
We use a velocity dispersion for the "halo" component of 
' H I Jialo = o-Hi,thin{hz,halo/hz,thin) (e.g., Kulkarni & Fich 
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1985; Combes & Becquaert 1997). Here ajji is the veloc- 
ity dispersion in the thin disk (7 km s -1 in all models) and 
h z .haio and h z ,thin are the scale height of the halo and the 
mean scale height of the flared thin disk, respectively. 

In Figures 3 and 4 we show the results for a model 
where the halo component has the following parame- 
ters: a FWHM of 3.3 kpc, a central volume density 16 
times lower than that of the disk material, a velocity dis- 
persion of 35 km s -1 , and a velocity that declines lin- 
early to 50% of the in-disk value at \z\ —2 kpc, and re- 
mains constant thereafter. This corresponds to a falloff of 
dV{r)/dz <25 km s _1 kpc -1 . While the particular combi- 
nation of model parameters we have used is not necessarily 
unique, we see from these figures that such a model rea- 
sonably reproduces the spatial and velocity distribution of 
Hi as observed in UGC 7321. The total intensity distri- 
bution derived from this model also matches well with the 
data (Figures 5 & 6). 

Given the success of the lagging halo model at repro- 
ducing the Hi structure of UGC 7321, we also explored 
whether the emission could be reproduced with a smooth 
disk and a lagging halo alone, without any flaring (the 
"smooth-l- lagging halo" model). The channel maps and 
vertical P-V slices predicted by such a model provide fair 
agreement with the data (Figures 3 & 4), although no- 
ticeable discrepancies begin to occur in the —98.2 km s -1 
channel map in Figure 4 and the 3?4W P-V slice in Fig- 
ure 3. Moreover, such a model clearly is seen to break 
down in Figures 5 & 6. One-dimensional vertical inten- 
sity slices from this model show more pronounced non- 
Gaussian wings than does the warp+flare+lagging halo 
model. And the total intensity image created from this 
model is seen to significantly underestimate the thickness 
of the disk at large projected radii. We conclude that 
UGC 7321 cannot be modelled without including a flaring 
of the gas layer. 

Lastly, we have also tested a model that includes warp- 
ing and flaring, but where the halo corotates with the disk 
(the "warp+flare+corotating halo" model). As seen in 
Figures 3 & 4, in several channels this model seems to pre- 
dict too much high-latitude emission at small projected 
radii (compare the V =56.8, —46.5, and —98.2 km s -1 
channel maps in Figure 4). However, only the outermost 
(~ 2d) contours are affected, and it is only at these lev- 
els that the differences between the lagging and corotating 
models become readily apparent. To complicate matters 
further, in some instances, vertically extended gas is seen 
at locations more consistent with the corotating model 
than the lagging halo model (for example, the vertical P-V 
cut at V.9 W in Figure 3). This could reflect noncircular 
motions of the emitting gas, or may indicate that the ro- 
tational lag of the halo has a dependence on r different 
from that in our models. At present, we are unable to 
distinguish between these possibilities. 

3.3.3. Tests of the Models in the Presence of Noise 

To better gauge how well we can distinguish between 
a lagging versus a corotating halo model given the noise 
in our observations, we have added to our models Gaus- 
sian noise identical to that in the real data and then re- 
examined the resulting channel maps and vertical P-V 
slices. We repeated this exercise several times in order 
to gauge how much faint features will appear to "vary" 



in the presence of random noise. We have also performed 
a similar experiment on the pure warp+flare model. Ex- 
amples of channel maps and vertical P-V slices from these 
"noisy" models are shown in Figures 7 & 8. 

From Figures 7 & 8 we can clearly see that noise alone is 
insufficient to emulate the types of high-latitude features 
we see in UGC 7321, as the warp+flare model shows es- 
sentially no evidence for such features. We also find that 
some differences between the lagging and corotating mod- 
els do persist even in the presence of noise. For example, 
in the channel maps over the velocity range 20-50 km s , 
the lagging model consistently shows high-latitude mate- 
rial "trailing" to large projected radius, behind the bulk of 
the low-latitude emission. Such features are consistently 
seen in the data, but are not reproduced by the corotat- 
ing halo model. Overall, at intermediate radii, the lagging 
halo model predicts more high-z emission near the sys- 
temic velocity than does the corotating halo model, and 
this also seems borne out by the data at most disk loca- 
tions where high-z Hi is observed (e.g., the 51", 78", and 
105"E P-V slices in Figure 8). 

Based on the above arguments, we suggest a tentative 
preference for the lagging halo model for UGC 7321 over 
one where the halo corotates with the disk, with the caveat 
that a lagging halo model would seem to require the addi- 
tion of some non-circular motions and/or a possible addi- 
tional change in V(z) as a function of r. However, given 
the limitations of present data and models, a situation 
where the bulk of the high-latitude gas corotates with the 
disk cannot yet be strongly ruled out. 

A further outcome of the exercise of introducing noise 
into our models is to illustrate that a smooth Hi halo 
can appear to break up into semi-discrete features in the 
presence of noise. Thus while the high-latitude gas in 
UGC 7321 may in fact have an intrinsically clumped or fil- 
amentary nature, this becomes difficult to distinguish from 
a relatively smooth emission distribution in our present 
data. 

4. DISCUSSION: THE ORIGIN OF THE VERTICALLY 
EXTENDED HI EMISSION 

We have shown in the previous sections that the LSB 
spiral UGC 7321 contains vertically extended Hi gas 
whose z-height extends beyond that predicted for a sim- 
ple, cool Hi layer, even after accounting for the effects of 
warping and flaring. Moreover, based on our best-fitting 
3-D models, we have suggested that at least a fraction of 
this high-latitude material may lag in rotation relative to 
the material in the midplanc. The caveats are that such 
a model would seem to require the addition of some non- 
circular motions and/or a possible additional change in 
V(z) as a function of r; also, given the limitations of our 
data and models, a situation where the bulk of the high- 
latitude gas corotates with the disk cannot yet be defini- 
tively ruled out. The emergence of a corotating halo model 
would have interesting implications for the origin of the 
high-latitude gas in UGC 7321, since both ballistic models 
(i.e., those associated with galactic fountain flows) and hy- 
drostatic models for galactic halos predict rotational lags. 
In the case of the hydrostatic models, cylindrical rotation 
can be achieved only if there exist additional pressure gra- 
dients or magnetic tension in the halo (Benjamin 2002). 
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4.1. Inferred Properties of the HI Halo 

Independent of the exact kinematic structure of the 
high-latitude gas in UGC 7321, we can use our models 
to compute an estimate of the total mass contained in the 
high- latitude Hi component in UGC 7321. For either our 
lagging or corotating halo model, the total mass of Hi 
in the halo component is ~ 1.3 x 10 8 Mq (roughly 12% 
of the total Hi content of UGC 7321). Of this material, 
- 1 x 10 7 M Q is present at \z\ > 30" (roughly 1% of the Hi 
content of UGC 7321). This latter value is approximately 
fifteen times less than the fractional Hi content inferred 
to be at comparable z-heights in NGC 891 by Swaters 
et al. (1997). We now briefly consider possibilities for 
the origin of this (comparatively weak) halo component in 
UGC 7321. 

4.2. Supernova Heating/ Galactic Fountain 

For those nearby galaxies whose gaseous halos have been 
studied to date, there appear to be well-established links 
between the properties of the multi-phase halos and the 
current star formation rates of the galaxies, implying the 
two are linked (e.g., Dahlcm 1997 and references therein). 
This suggests that at least some portion of the halos may 
be formed via internal processes. 

The type of thickened, rotationally lagging H I layer that 
we infer for UGC 7321 in Section 3.3 is qualitatively similar 
to the expectation for material originating via a galactic 
fountain, as described by Bregman (1980). In this sce- 
nario, gas heated by supernovae is predicted to move both 
toward larger \z\ and larger r before cooling, recondensing, 
and returning near the point of origin via ballistic orbits. 
While in the halo, the gas should lag the disk rotation as a 
consequence of conservation of angular momentum. Using 
qualitative arguments, Swaters et al. (1997) and Schaap 
et al. (2000) have suggested scenarios of this type may 
account for the origin of the Hi halos of NGC 891 and 
NGC 2403, respectively. However, it is unclear whether a 
detectable H I halo could be formed in this manner in a 
galaxy like UGC 7321, that exhibits such modest levels of 
current star formation. 

Based on the IRAS far-infrared 60 and lOOyum fluxes 
for UGC 7321 (0.344 Jy±14% and 0.964 Jy±19%, 
respectively) 6 and the formulae given by Condon (1992), 
we estimate the current star formation rate in UGC 7321 
for stars with M > 5 Mq is only ^0.006 Mq per year. 
This is ^200 times smaller than the value similarly derived 
for NGC 891. On the other hand, UGC 7321 is a smaller, 
lower mass spiral than others where high-latitude H I has 
been reported, and likely also has a more tenuous inter- 
stellar medium. These factors may increase the efficiency 
with which supernovae can deposit gas at high latitudes. 

An independent estimate of the star formation rate in 
UGC 7321 can be derived from the Ha data of Matthews 
et al. (1999). Although the Ha+[Nn] image presented by 
these authors is not flux calibrated, using the known R- 
band luminosity of the galaxy, the mean Ha / [N n] line ra- 
tio from Goad & Roberts (1981), and the relative through- 
puts of the broad-band continuum (R) and narrow-band 
(Ha+[Nn]) filters used by Matthews et al., we estimate 
an Ha luminosity for UGC 7321 from their data of Ln a ~ 
1 x 10 40 erg s _1 . Such a value is typical of quiescent, late- 

6 Valucs taken from the NED database. 



type spirals (e.g., Thilker et al. 2002). Using the formula 
given by Condon (1992), this corresponds to a massive star 
formation rate of ^0.02 Mq per year-a few times higher 
than estimated above from the FIR emission, but roughly 
consistent to within the uncertainties of both estimates. 

In order to test whether it is feasible that supernovae 
could supply the energy required to maintain an observable 
Hi halo in UGC 7321, we follow an approach similar to 
Kulkarni & Fich (1985). We neglect possible partial sup- 
port for the gas that may come from magnetic fields, cos- 
mic rays, and radiation, since based on the extremely weak 
radio continuum flux from this galaxy (Uson & Matthews 
2003), these terms are likely to be quite small. In this 
case, assuming hydrostatic equilibrium, the kinetic energy 
of the halo is simply given by E kMlo = ^M halo a 2 HI halo , 
where Mhaio is the mass of the gaseous halo material (H I 
plus He) and a hi, halo is the velocity dispersion of the halo 
(35 km s _1 for our model). Using the halo mass quoted in 
Section 4.1, applying a factor of 1.34 to correct for He, and 
assuming a radius for the galaxy of 11.8 kpc, we derive the 
kinetic energy per unit area of Ek,haio/A = 5.1 x 10 s erg 
cm~ 2 for the gaseous halo. 

We assume that the timescale for intercloud collisions in 
the halo is likely to be longer than their free-fall time back 
to the disk, and therefore use the latter timescale to esti- 
mate the rate at which energy must be continuously sup- 
plied to maintain the halo. Assuming a downward velocity 
of 35 km s" 1 from a height of 1.6 kpc, the approximate 
free-fall time is tff ~46 million years. This then requires 
an energy input rate of E kha i /{Attt) w 3.6 x 10~ 7 erg 

— 1—2 

s cm . 

From Cappellaro et al (1993), we find that mean su- 
pernova rates for typical Scd/Sd galaxies are 2.43±0.70 
per 100 years, per 10 10 i Q , of which ^24% are typically 
type la and ^76% are core collapse supernovae. Since the 
B-band luminosity of UGC 7321 is l.lxlO 9 i (Uson & 
Matthews 2003), this yields an estimated total supernova 
rate of ~ 3 x 10~ 3 year -1 . This number is roughly a fac- 
tor of two higher than we might estimate from the massive 
star formation rate derived from the FIR fluxes above, but 
given the large uncertainty in the FIR measurements, we 
adopt this as a more "optimistic" estimate. 

Fiducial kinetic energies for type la and core collapse 
supernovae are 5 x 10 50 ergs and 10 51 ergs, respectively 
(Abbott 1982). Of this energy, approximately 1% and 
3%, respectively can be converted into bulk motions of gas 
clouds (Spitzer 1978). This predicts a kinetic energy per 
unit area supplied by supernovae of E$n/A ~ 5.0 x 10~ 7 
erg s" 1 cm -2 — a value comparable to Ek,haio/(Atff) esti- 
mated above, and therefore consistent with the possibility 
that supernova energies may be sufficient to maintain the 
high-latitude Hi observed in UGC 7321. 

We can provide an independent check on our energy in- 
jection estimate using the figures and relations supplied by 
Collins et al. (2002). Using their Figure 4, we find that a 
maximum observed gas height of 2.4 kpc implies a ratio of 
"kick" velocity V% to circular velocity V c in the midplanc 
of ~0.8. Choosing some intermediate galactocentric radius 
(say 4 kpc), and taking V(r) «100 km s" 1 for UGC 7321, 
this implies 14=80 km s" 1 . From Figure 6 of Collins et al., 
the implied "cycling frequency" is f cyc i e = 1.4 x 10 _8 yr _1 . 
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Inserting this into their Equation 5 leads to an esti- 
mate of the halo mass flux of Mh ~ 2.3 M Q yr" 1 . Fi- 
nally, from the relation for the energy input requirement: 
E = (1.1 x lO^ergs-^MhV^oo (where the halo mass flux 
is in units of M Q yr _1 and the kick velocity is in units 
of 100 km s _1 ) we arrive at E = 1.6 x 10 39 erg s _1 , or 
E/A = 3.8 x 10~ 7 erg s _1 cm -2 — very close to our esti- 
mate above. 

Using the "kick" velocity estimated from the Collins 
et al. models, we can make another interesting compar- 
ison, namely how that value compares with the escape 
speed from the galaxy, v e . Assuming a spherical poten- 
tial, the escape speed from a galaxy is given by equa- 
tion 2-292 of Binney & Trcmainc (1987). At our fidu- 
cial radius of r =4 kpc, we thus derive an escape speed 
v e ssl41 km s _1 — significantly higher than the kick veloc- 
ity of 80 km s _1 , implying UGC 7321 would be able to 
retain most of its fountain gas. 

Of course the above estimates are highly idealized, and 
more sophisticated computations are clearly needed to test 
more accurately how efficiently supernovae (and stellar 
winds) can be in supplying gas to the halo, and whether 
the observed spatial and velocity structures of the high- 
latitude material quantitatively matches the predictions of 
galactic fountain-type models. Very deep Ha observations 
to search for an extraplanar component of ionized gas in 
UGC 7321 would also be of interest in further testing the 
link between star formation and high-latitude Hi in this 
galaxy. In the mean time, we also consider the feasibility 
of other possibilities for the origin of the halo material in 
UGC 7321. 

4.3. Alternative Possibilities 

In spite of the well-established links between current star 
formation and the existence of gaseous halos in brighter 
galaxies, it is possible that star formation simply acts to 
augment the halo and lead to the partial ionization of pre- 
existing halo gas, thereby making these halos more exten- 
sive and more readily observable via a variety of tracers. 
This may be expected to occur if some gas were deposited 
in the halo either during an earlier phase in the galaxy's 
evolution, or via ongoing accretion. We now consider these 
two possibilities in turn for the source of the high-latitude 
Hi in UGC 7321. 

4.3.1. Galactic Superwinds 

Steidel et al. (2002) recently reported that for a sample 
of five intermediate redshift Mgll absorption line galax- 
ies, the models needed to explain the kinematic properties 
of the absorbers require either thick rotating gas layers 
and/or rotational velocities that decline with height above 
the plane. One galaxy in their sample, Gl 1222+228, 
is particularly interesting in light of our findings for 
UGC 7321, as it appears to show some intriguing analo- 
gies. 

Like UGC 7321, Gl 1222+228 (redshift 0.55) has the 
morphology of a pure disk system of relatively low surface 
brightness, in spite of its edge-on geometry. Its absolute 
blue magnitude (—18.09) and its peak rotational veloc- 
ity (100 km s _1 ) both are also very similar to UGC 7321. 
From Stcidcl et al.'s Figure 4, we estimate the linear di- 
ameter of Gl 1222+228 to be <~12.1/i _1 kpc — again com- 



parable to UGC 7321, particularly if the faint outskirts of 
the disk are invisible owing to their LSB nature. For this 
system, Steidel et al. successfully reproduced the Mgll 
absorption profile for this galaxy using a model with an 
"effective height" for the gas layer (as measured from the 
midplanc) of <5 kpc and a velocity structure for the high- 
z-height material that falls to zero velocity within a few 
kpc above the disk plane. This model is interestingly sim- 
ilar to our lagging halo model for UGC 7321. 

Steidel et al. (2002) note the general problem of ex- 
plaining Mgll halos of moderate redshift galaxies from 
galactic fountains, since most Mgll absorption systems 
are not currently highly active star-formers. Motivated by 
results from the strong far-UV absorption lines observed 
in Lyman break galaxies, they suggest as an alternative, 
large-scale galactic winds from a much earlier epoch of star 
formation (at rcdshifts >3). In this scenario, the shock- 
heated and metal-enriched supcrwind material would cool 
slowly to form the gaseous disk and halo observed in nor- 
mal disk galaxies at low or moderate redshifts. Could this 
be the origin of the high-latitude Hi in UGC 7321? 

Although the galactic superwind scenario suggested by 
Stcidcl ct al. (2002) has not yet been quantified in detail, 
based on dynamical arguments it is not clear whether ma- 
terial deposited in the halo at high redshift would still be 
observable in local galaxies. The metal-line absorption fea- 
tures in the halos of galaxies probed by QSO absorption 
lines, including Mgll, appear to arise in discrete clouds 
(e.g., Petitjean & Bergeron 1990), and Chen, Lanzetta, 
& Webb (2001a) have argued that the timescale for colli- 
sions and subsequent orbital decay of these clouds should 
be considerably less than a Hubble time. Therefore it is 
uncertain whether any remnant of gaseous halos formed by 
superwinds at high rcdshifts could persist to the present 
epoch. 

4.3.2. Inf ailing Gas 

If the high-latitude Hi material in UGC 7321 is not di- 
rectly tied to current star formation, and is not a remnant 
from an early galactic superwind, an alternative is that at 
least some fraction of this gas has been acquired through 
the accretion of intergalactic clouds or small companions 
(e.g., Oort 1970; van der Hulst & Sancisi 1988). This 
could have occurred in the form of a one-time accretion 
of a dwarf satellite or the "drizzle" of small parcels of H I 
onto the galaxy. 

UGC 7321 appears to be a very isolated galaxy, and 
Uson & Matthews (2003) put a lOcr upper limit on the H I 
masses of any gas-rich companions or clumps within 12 
arcminutes (36 kpc) of UGC 7321 of < 2.2 x 10 6 M Q . In 
addition, the mass of any possible infalling clumps would 
be limited by the constraint that they cannot excessively 
heat and disrupt the disk (e.g., Toth & Ostriker 1992). 

We can obtain a rough estimate of the maximum mass 
of a satellite, M s , that could have been recently accreted 
by UGC 7321 without destroying its thin stellar disk using 
the relation given by Lacey (1991): 

Ms <ol*M D V~ 2 ax . 

Here, cr z .* is the stellar velocity dispersion, Mo is the 
disk mass, and V max is the peak disk rotational veloc- 
ity. Taking Mb to be equal to the total estimated mass 
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of gas and stars in the disk, V ma x = 105 km s 1 (Uson & 
Matthews 2003), and a z , *=20 km s" 1 (Matthews 2000), 
we find M s ^7.5 x 10 7 Mq. This is comparable to our 
inferred halo mass for UGC 7321; however, any gas-rich 
dwarf galaxy is likely to have a dark matter content equal 
to at least a few times its visible mass (e.g., Cote, Carig- 
nan, & Freeman 2000). Thus it appears possible, although 
unlikely that the high-latitude material in UGC 7321 orig- 
inated from a single accretion event. 

Still another possibility is that UGC 7321 has been ac- 
creting small Mhi < 2 x 10 6 M Q parcels of gas throughout 
much of its lifetime. Indeed, a variety of lines of evidence 
suggest that slow, continual gas accretion (of order a few 
Mq per year) may be a common occurrence for disk galax- 
ies (e.g., Wakker et al. 1999; Binney 2000; see also Oort 
1970). In the case of UGC 7321, the misaligned angu- 
lar momentum thus acquired (e.g., Binney & Jiang 1999; 
Dcbattista & Sellwood 1999) could also offer an explana- 
tion for the origin of the warp in this isolated system (see 
also the discussion in Uson & Matthews 2003). In order 
to maintain the observed thick H I layer, we estimate that 
a mass infall rate of Jj3.8M Q per year would be required 

based on the relation E — jMV 2 , where we have taken 

E = 1.5 x 10 39 erg s" 1 and assumed V «35 km s _1 . This 
is comparable to gas infall rates estimated for the Milky 
Way (Binney 2000). However, in order for infalling mate- 
rial to share the rotation of the disk, this requires that the 
clumps must become "dragged" along with some sort of 
pre-existing, rotating medium, and the existence of such 
a medium is unclear. This objection could partially be 
alleviated if infalling clumps acted to trigger further in- 
stabilities in the disk and thus raised additional rotating 
gas to high z. 

5. SUMMARY 

Using sensitive VLA Hi observations in combination 
with 3-D modelling, we have found evidence for high- 
latitude (\z\ <2.4 kpc) Hi emission in the edge-on, LSB 
spiral galaxy UGC 7321. We are unable to reproduce the 
spatial and kinematic structures of this vertically extended 
material solely by projection effects from warping and flar- 
ing of the gas layer. However, we are successfully able to 
model the Hi distribution in UGC 7321 with a combina- 



tion of a thin, cool H I disk with mild warping and flaring, 
together with a thicker "halo" component. 

We find tentative evidence that the vertically ex- 
tended gas exhibits a lag in rotation velocity by of 
<;25 km s _1 kpc -1 relative to the cooler Hi material in 
the midplane. However, given the possible presence of ad- 
ditional non-circular motions and the sensitivity limits of 
our data, it is impossible to completely rule out models 
where the extended gas largely corotates with the disk. 

We estimate a FWHM for the cool disk component to be 
^260 pc at the disk center, flaring to a FWHM of 1.6 kpc 
at the last observed point; we estimate the FWHM thick- 
ness of the halo to be ~3.3 kpc. Based on our best models, 
the Hi "halo" comprises roughly 12% of the total Hi con- 
tent of UGC 7321. Although the best-fitting model we 
have derived for UGC 7321 is not necessarily unique, af- 
ter exploring a large range of parameter space, we find 
our model to reproduce the data significantly better than 
other physically motivated models, including those with a 
higher-amplitude flare, or those with a vertically extended 
halo but no flaring. 

Our detection of vertically extended Hi emission (be- 
yond that attributable to a warp or flare) in UGC 7321 
represents the first time that evidence of such an ISM com- 
ponent has been reported in a nearby galaxy with such a 
low current star formation rate. Our rough calculations 
imply that energy input from supernovae may be sufficient 
to maintain the observed high-latitude Hi in UGC 7321. 
However, slow, sustained gaseous infall (in the form of 
clumps with Hi masses ^10 6 Mq) remains another vi- 
able explanation, and could also explain the existence of 
the warp in this isolated galaxy. Our findings suggest that 
H I halos may be common features of spiral disk galaxies in 
the local universe, and that their existence is not necessar- 
ily dependent on vigorous levels of current star formation. 

We are grateful to Jay Lockman, Bob Benjamin, Mark 
Rcid, and Hsiao- Wen Chen for valuable suggestions and 
discussions, and to an anonymous referee whose comments 
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from a Clay Fellowship from the Harvard-Smithsonian 
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Fig. 1. — Hi total intensity map for UGC 7321 overlayed on an ij-band image from Matthews et al. 1999. Contour levels are 
2,2.8,4,. ..44,64,88x(7.7 X 10 19 ) atoms cm" 2 . The FWHM beam size is ~16". 
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Fig. 2. — Hi intensity profiles extracted perpendicular to the midplanc of UGC 7321 at several locations along the disk using the data 
shown in Figure 1 (thin solid lines). Single Gaussian fits to the data are overplotted as dotted lines. FWHM of the fitted Gaussians range 
from 22'.'4-24'.'4. For comparison, the thick solid line on the center panel shows a Gaussian with FWHM equal to the synthesized beamwidth 
(16"). At all locations, extended wings are visible in the data in excess of the Gaussian models. 
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Fig. 3. — Position-velocity profiles extracted parallel to the minor axis of UGC 7321 at several disk locations (right column) and at the 
corresponding locations in various galaxy models discussed in Section 3.3. Axes are velocity (in km s -1 with respect to the systemic velocity) 
and the distance from the disk midplane (in arcseconds). Contour levels are — 0.2,0.2, 0.4,0. 8,.... 12. 8x (0.39 mjy beam -1 ). The lowest contour 
level is ~ 2a. The profiles shown in each row were extracted at 1.'5 intervals: 2.'6 and l.'l east of the minor axis; 0.'4, 1.'9, and 3.'4 west of the 
minor axis. 




Fig. 4. — Representative Hi channel maps from our UGC 7321 VLA data (righthand column) and several of the galaxy models discussed in 
Section 3.3. The upper two panels in each column are from the receding side of the disk, while the lower four panels are from the approaching 
side. Contour levels are -0. 4, -0.2, 0.2, 0.4. ..6. 4, lOx (0.39 mjy beam -1 ) (spatial resolution ~ 16" and velocity resolution ~5.2 km s~ 1 ), with 
the lowest contour approximately 2a. 
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High-Latitude Hi in UGC 7321 
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Fig. 5. — Same as Figure 2, but for the "warp+flare" , "warp+flare+lagging halo", "smooth+lagging halo", and "warp+flare+corotating 
halo" models. 
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Fig. 6. — Hi total intensity maps for the various galaxy models discussed in Section 3.3 (lower panels) and for our data (upper panel). 
Contour levels are 1.4,2,2.8,4,...44,64,88x (7.7 X 10 19 ) atoms cm~ 2 (spatial resolution^ 16" and velocity resolution ~5.2 km s _1 ). Only points 
in each channel map with flux > — 1.5a were included when computing in these images. 
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High-Latitude Hi in UGC 7321 
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Fig. 7. — Selected channel maps from three of our models with Gaussian noise added (see Section 3.3.3). The corresponding channel maps 
from the real data are also shown. Contour levels are as in Figure 4. 
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Fig. 8. — Position-velocity profiles extracted parallel to the minor axis at several disk locations for three of our models, with Gaussian noise 
added. The real data at the corresponding positions are also shown. The profiles in each row were extracted at 27" intervals between 159" 
east and 24" east of the minor axis. Contour levels are as in Figure 3. 



